Estimating the Effective Reproduction Number of COVID-19 in United Kingdom
Louis Rossouw
2021-03-30 22:45:02 BST
1 Introduction
This paper contains estimates for the effective reproduction number \(R_{t,m}\) over time \(t\) in various nations and other regions \(m\) of United Kingdom. This is done using the methodology as described in [1]. These have been implemented in R using EpiEstim package [2] which is what is used here. The methodology and assumptions are described in more detail here.
This paper and it’s results should be updated roughly daily and is available online.
As this paper is updated over time this section will summarise significant changes. The code producing this paper is tracked using Git. The Git commit hash for this project at the time of generating this paper was 468c0059dfbff1fe11d7ecdd2ca657932bde77b8.
2 Data
Data is obtained [3]. This contains the daily cases, hospital admissions and deaths for United Kingdom by various geographies. Here the data is accessed by specimen date, admission date and date of death.
Once history is built up an allowance for rate reported cases, admissions and deaths will be made. For now the data is cut-off a week prior to the last case date in the data. All data prior to 1 March 2020 are removed.
3 Methodology
The methodology is described in detail here.
4 Results by Nation
4.1 Cases
Below a 7-day moving average of daily case count on a log scale by nation is plotted:
Daily Cases by Nation (7-day moving average)
4.2 Admissions
Below we plot cumulative hospital admissions on a log scale.
Daily Admissions by Nation (7-day moving average)
4.3 Deaths
Below a 7-day moving average of daily deaths by nation on a log scale is plotted:
Daily Deaths by Nation (7-day moving average)
4.4 Current \(R_{t,m}\) estimates by Nation
Below current (last weekly) \(R_{t,m}\) estimates are tabulated.
| Nation | Estimate Type | Count (Last Week) | Week Ending | R - Lower CI | R - Mean | R - Uppper CI |
|---|---|---|---|---|---|---|
| England | cases | 31,610 | 2021-03-23 | 1.0 | 1.0 | 1.0 |
| England | admissions | 2,132 | 2021-03-23 | 0.8 | 0.8 | 0.9 |
| England | deaths | 374 | 2021-03-23 | 0.6 | 0.7 | 0.8 |
| Northern Ireland | cases | 1,034 | 2021-03-23 | 0.9 | 0.9 | 1.0 |
| Northern Ireland | admissions | 58 | 2021-03-23 | 0.6 | 0.8 | 1.0 |
| Northern Ireland | deaths | 9 | 2021-03-23 | 0.5 | 1.1 | 1.9 |
| Scotland | cases | 3,901 | 2021-03-23 | 1.0 | 1.0 | 1.0 |
| Scotland | admissions | 186 | 2021-03-23 | 0.8 | 0.9 | 1.0 |
| Scotland | deaths | 28 | 2021-03-23 | 0.5 | 0.7 | 1.0 |
| Wales | cases | 1,231 | 2021-03-23 | 0.9 | 0.9 | 1.0 |
| Wales | admissions | 332 | 2021-03-23 | 1.0 | 1.1 | 1.3 |
| Wales | deaths | 21 | 2021-03-23 | 0.4 | 0.7 | 1.0 |
Estimated Effective Reproduction Number by Nation
4.5 Maps of Effective Reproduction Number
Below we plot the current effective reproduction number on maps with boundaries provided by [4].
4.5.1 Cases
4.5.2 Admissions
4.5.3 Deaths
4.6 Graphs over Time
Below we plot results for each nation We filter out weeks where the upper end of confidence interval for \(R_{t,m}\) exceeds 4.
4.6.1 England
4.6.2 Northern Ireland
4.6.3 Scotland
4.6.4 Wales
5 Results by Region
5.1 Cases
Below we daily case count is plotted on a log scale by region:
Daily Cases by Region (7-day moving average)
5.2 Deaths
Below a 7-day moving average of daily deaths by region is plotted on a log scale:
Daily Deaths by Region (7-day moving average)
5.3 Current \(R_{t,m}\) estimates by Region
Below current (last weekly) \(R_{t,m}\) estimates are tabulated.
| Region | Estimate Type | Count (Last Week) | Week Ending | R - Lower CI | R - Mean | R - Uppper CI |
|---|---|---|---|---|---|---|
| Northern Ireland | cases | 1,034 | 2021-03-23 | 0.9 | 0.9 | 1.0 |
| Northern Ireland | deaths | 9 | 2021-03-23 | 0.5 | 1.1 | 1.9 |
| Scotland | cases | 3,901 | 2021-03-23 | 1.0 | 1.0 | 1.0 |
| Scotland | deaths | 28 | 2021-03-23 | 0.5 | 0.7 | 1.0 |
| Wales | cases | 1,231 | 2021-03-23 | 0.9 | 0.9 | 1.0 |
| Wales | deaths | 21 | 2021-03-23 | 0.4 | 0.7 | 1.0 |
| East Midlands | cases | 3,673 | 2021-03-23 | 1.0 | 1.0 | 1.0 |
| East Midlands | deaths | 58 | 2021-03-23 | 0.6 | 0.8 | 1.0 |
| East of England | cases | 2,751 | 2021-03-23 | 1.0 | 1.0 | 1.0 |
| East of England | deaths | 31 | 2021-03-23 | 0.4 | 0.6 | 0.8 |
| London | cases | 3,258 | 2021-03-23 | 0.9 | 1.0 | 1.0 |
| London | deaths | 20 | 2021-03-23 | 0.4 | 0.6 | 0.9 |
| North East | cases | 1,779 | 2021-03-23 | 0.9 | 0.9 | 1.0 |
| North East | deaths | 14 | 2021-03-23 | 0.4 | 0.7 | 1.2 |
| North West | cases | 5,041 | 2021-03-23 | 0.9 | 1.0 | 1.0 |
| North West | deaths | 70 | 2021-03-23 | 0.6 | 0.8 | 1.0 |
| South East | cases | 3,155 | 2021-03-23 | 0.9 | 1.0 | 1.0 |
| South East | deaths | 42 | 2021-03-23 | 0.4 | 0.6 | 0.8 |
| South West | cases | 1,659 | 2021-03-23 | 1.0 | 1.0 | 1.1 |
| South West | deaths | 20 | 2021-03-23 | 0.4 | 0.7 | 1.1 |
| West Midlands | cases | 3,791 | 2021-03-23 | 0.9 | 1.0 | 1.0 |
| West Midlands | deaths | 57 | 2021-03-23 | 0.6 | 0.8 | 1.0 |
| Yorkshire and The Humber | cases | 6,297 | 2021-03-23 | 1.0 | 1.0 | 1.1 |
| Yorkshire and The Humber | deaths | 53 | 2021-03-23 | 0.5 | 0.7 | 0.9 |
Estimated Effective Reproduction Number by Region
5.4 Maps of Effective Reproduction Number
Below we plot the current effective reproduction number on maps with boundaries provided by [5].
5.4.1 Cases
5.4.2 Deaths
5.5 Graphs over Time
Below we plot results for each nation.
5.5.1 Northern Ireland
5.5.2 Scotland
5.5.3 Wales
5.5.4 East Midlands
5.5.5 East of England
5.5.6 London
5.5.7 North East
5.5.8 North West
5.5.9 South East
5.5.10 South West
5.5.11 West Midlands
5.5.12 Yorkshire and The Humber
6 Results by NHS Region
6.1 Current \(R_{t,m}\) estimates by Region
Below current (last weekly) \(R_{t,m}\) estimates are tabulated.
| Region | Estimate Type | Count (Last Week) | Week Ending | R - Lower CI | R - Mean | R - Uppper CI |
|---|---|---|---|---|---|---|
| Northern Ireland | admissions | 58 | 2021-03-23 | 0.6 | 0.8 | 1.0 |
| Scotland | admissions | 186 | 2021-03-23 | 0.8 | 0.9 | 1.0 |
| Wales | admissions | 332 | 2021-03-23 | 1.0 | 1.1 | 1.3 |
| East of England | admissions | 207 | 2021-03-23 | 0.7 | 0.8 | 0.9 |
| London | admissions | 225 | 2021-03-23 | 0.7 | 0.8 | 0.9 |
| North West | admissions | 377 | 2021-03-23 | 0.8 | 0.9 | 1.0 |
| South East | admissions | 254 | 2021-03-23 | 0.7 | 0.8 | 0.9 |
| South West | admissions | 128 | 2021-03-23 | 0.6 | 0.8 | 0.9 |
| Midlands | admissions | 465 | 2021-03-23 | 0.7 | 0.8 | 0.9 |
| North East and Yorkshire | admissions | 476 | 2021-03-23 | 0.8 | 0.8 | 0.9 |
Estimated Effective Reproduction Number by Region
6.2 Maps of Effective Reproduction Number
Below we plot the current effective reproduction number on maps with boundaries provided by [6].
6.3 Graphs over Time
Below we plot results for each nation We filter out weeks where the upper end of confidence interval for \(R_{t,m}\) exceeds 4.
6.3.1 East of England
6.3.2 London
6.3.3 Midlands
6.3.4 North East and Yorkshire
6.3.5 North West
6.3.6 South East
6.3.7 South West
7 Results by Upper Tier Local Authority
7.1 Highest \(R_{t,m}\) as estimated using cases
Below we plot Upper Tier Local Authorities with the highest reproduction numbers (providing they had at least 100 cases in last 7 days):
## Selecting by Rt_ui_95| Upper Tier Local Authority | Estimate Type | Count (Last Week) | Week Ending | R - Lower CI | R - Mean | R - Uppper CI |
|---|---|---|---|---|---|---|
| Moray | cases | 68 | 2021-03-23 | 1.9 | 2.6 | 3.3 |
| Barking and Dagenham | cases | 141 | 2021-03-23 | 1.3 | 1.6 | 1.9 |
| Bexley | cases | 96 | 2021-03-23 | 1.2 | 1.5 | 1.8 |
| Derry City and Strabane | cases | 132 | 2021-03-23 | 1.1 | 1.4 | 1.7 |
| Greenwich | cases | 140 | 2021-03-23 | 1.1 | 1.3 | 1.5 |
| Hackney and City of London | cases | 64 | 2021-03-23 | 1.0 | 1.3 | 1.6 |
| Dundee City | cases | 151 | 2021-03-23 | 1.1 | 1.3 | 1.5 |
| West Dunbartonshire | cases | 72 | 2021-03-23 | 1.0 | 1.3 | 1.6 |
| Ards and North Down | cases | 125 | 2021-03-23 | 1.0 | 1.3 | 1.5 |
| Aberdeen City | cases | 98 | 2021-03-23 | 1.0 | 1.3 | 1.5 |
| Plymouth | cases | 112 | 2021-03-23 | 1.0 | 1.2 | 1.5 |
| Windsor and Maidenhead | cases | 89 | 2021-03-23 | 1.0 | 1.2 | 1.5 |
| Central Bedfordshire | cases | 157 | 2021-03-23 | 1.0 | 1.2 | 1.4 |
| Aberdeenshire | cases | 100 | 2021-03-23 | 1.0 | 1.2 | 1.4 |
| Blackpool | cases | 67 | 2021-03-23 | 0.9 | 1.2 | 1.5 |
| Swansea | cases | 154 | 2021-03-23 | 1.0 | 1.2 | 1.4 |
| Enfield | cases | 114 | 2021-03-23 | 1.0 | 1.2 | 1.4 |
| Trafford | cases | 119 | 2021-03-23 | 1.0 | 1.2 | 1.4 |
| East Dunbartonshire | cases | 63 | 2021-03-23 | 0.9 | 1.2 | 1.5 |
| Cornwall and Isles of Scilly | cases | 105 | 2021-03-23 | 0.9 | 1.2 | 1.4 |
| Telford and Wrekin | cases | 65 | 2021-03-23 | 0.9 | 1.1 | 1.4 |
| Hartlepool | cases | 94 | 2021-03-23 | 0.9 | 1.1 | 1.4 |
| Hammersmith and Fulham | cases | 68 | 2021-03-23 | 0.9 | 1.1 | 1.4 |
| Neath Port Talbot | cases | 72 | 2021-03-23 | 0.9 | 1.1 | 1.4 |
| Merthyr Tydfil | cases | 82 | 2021-03-23 | 0.9 | 1.1 | 1.4 |
7.2 Risk Quadrants
The plots below show weekly cases (or deaths) on the X-axis and the reproduction number on the Y-axis. By dividing this into 4 quadrants we can identify upper tier local authorities with high cases and high reproduction numbers, or high cases and low reproduction numbers etc.
Values where the reproduction number exceeds 3 are plotted at 3.
7.2.1 Cases
Risk Quadrants - Cases
7.2.2 Deaths
Risk Quadrants - Deaths
7.3 Map of Effective Reproduction Number (Cases)
Below we plot the current effective reproduction number estimated from case data on maps with boundaries provided by [7].
7.4 Map of Reproduction Number by Upper Tier Local Authority Over 60 days
Below the reproduction number by week for each Upper Tier Local Authority is animated over last 60 days:
8 Detailed Results
Detailed output are saved to a comma-separated value file. The file can be found here.
9 Discussion
Limitation of this method to estimate \(R_{t,m}\) are noted in [1]
- It’s sensitive to changes in transmissibility, changes in contact patterns, depletion of the susceptible population and control measures.
- It relies on an assumed generation interval assumptions.
- The size of the time window can affect the volatility of results.
- Results are time lagged with regards to true infection, more so in the case of the use of deaths.
- It’s sensitive to changes in case (or death) detection.
- The generation interval may change over time.
Further to the above the estimates are made under assumption that the cases and deaths are reported consistently over time. For cases this means that testing needs to be at similar levels and reported with similar lag. Should these change rapidly over an interval of a few weeks the above estimates of the effective reproduction numbers would be biased. For example a rapid expansion of testing over the last 3 weeks would results in overestimating recent effective reproduction numbers. Similarly any changes in reporting (over time and underreporting) of deaths would also bias estimates of the reproduction number estimated using deaths. It may well be that some catch-up in reported deaths is exaggerating the estimates for October.
Estimates for the reproduction number are plotted in time period in which the relevant measure is recorded. Though in reality the infections giving rise to those estimates would have occurred roughly between a week to 4 weeks earlier depending on whether it was cases or deaths. These figures have not been shifted back.
Despite these limitation we believe the ease of calculation of this method and the ability to use multiple sources makes it useful as a monitoring tool.
11 Digital boundaries
Office for National Statistics licensed under the Open Government Licence v.3.0
Contains OS data © Crown copyright and database right 2020
References
[1] A. Cori, N. M. Ferguson, C. Fraser, and S. Cauchemez, “A new framework and software to estimate time-varying reproduction numbers during epidemics,” American Journal of Epidemiology, vol. 178, no. 9, pp. 1505–1512, Sep. 2013, doi: 10.1093/aje/kwt133. [Online]. Available: https://doi.org/10.1093/aje/kwt133
[2] A. Cori, EpiEstim: A package to estimate time varying reproduction numbers from epidemic curves. 2013 [Online]. Available: https://CRAN.R-project.org/package=EpiEstim
[3] Office for National Statistics, “Official UK Coronavirus Dashboard,” 2020. [Online]. Available: https://coronavirus.data.gov.uk. [Accessed: 07-Nov-2020]
[4] Office for National Statistics, “Countries (December 2019) Boundaries UK BUC,” 09-Oct-2017. [Online]. Available: https://geoportal.statistics.gov.uk/search?collection=Dataset. [Accessed: 07-Nov-2020]
[5] Office for National Statistics, “NUTS Level 1 (January 2018) Ultra Generalised Clipped Boundaries in the United Kingdom,” 31-Jul-2017. [Online]. Available: https://geoportal.statistics.gov.uk/datasets/nuts-level-1-january-2018-ultra-generalised-clipped-boundaries-in-the-united-kingdom. [Accessed: 07-Nov-2020]
[6] Office for National Statistics, “NHS England Regions (April 2020) Boundaries EN BUC,” 13-May-2020. [Online]. Available: https://geoportal.statistics.gov.uk/search?collection=Dataset. [Accessed: 07-Nov-2020]
[7] Office for National Statistics, “Counties and Unitary Authorities (December 2019) Boundaries UK BUC,” 11-Mar-2020. [Online]. Available: https://geoportal.statistics.gov.uk/datasets/counties-and-unitary-authorities-december-2019-boundaries-uk-buc. [Accessed: 07-Nov-2020]